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^y. , Using the formalism of supersymmetric quantum mechanics, we obtain a large number of new 

^ ■ analytically solvable one-dimensional periodic potentials and study their properties. More specif- 

Q,' ically, the supersymmetric partners of the Lame potentials mala + 1) svP'(x,m) are computed for 
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C^ \ distinctly different from the original Lame potential, even though they both have the same energy 

Q^. band structure. We also derive and discuss the energy band edges of the associated Lame po- 
tentials pm ST?{x,m) -\- qm cn^(a;, TO)/dn [x^m), which constitute a much richer class of periodic 

^^ ' problems. Computation of their supersymmetric partners yields many additional new solvable and 

Cu ^ quasi exactly solvable periodic potentials. 
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1. Introduction: 

The energy spectrum of electrons on a lattice is of central importance in condensed matter physics. 
In particular, knowledge of the existence and locations of band edges and band gaps determines many 
physical properties. Unfortunately, even in one dimension, there are very few analytically solvable 
periodic potential problems in quantum mechanics. The aim of this paper is to extend the small 
currently known set of analytically solvable periodic potentials. 

For a potential with period L, one is seeking solutions of the Schrodinger equation subject to the 
Bloch condition 

V'(x) = e'^^ ij{x + L) , (1) 

where k denotes the crystal momentum. The spectrum shows energy bands whose edges correspond 
to kL = 0,7r, that is the wave functions at the band edges satisfy ip{x) = ib^(x + L). For periodic 
potentials, the band edge energies and wave functions are often called eigenvalues and eigenfunctions, 
and we will also use this terminology. The classic text book example which is used to demonstrate 
band structure is the Kronig-Penney model 

oo 

V{x) = J2 ^oS(^ - ^^) ■ 

n=—oo 

It should be noted that the band edges for the Kronig-Penney model can only be computed by solving 
a transcendental equation. Another well studied class of periodic potentials is 

V{x) = pm sv?{x, m) , p = a{a + 1) . (2) 

Here sn(x,m) is a Jacobi elliptic function of real elliptic modulus parameter m (0 < m < 1) with 
period AK{m). For simplicity, from now onward, we will not explicitly display the modulus parameter 
m as an argument of Jacobi elliptic functions |]l|. The elliptic function potentials of eq. (|2|) have a 
period L = 2K{m), and will be referred to as Lame potentials, since the corresponding Schrodinger 
equation is called Lame's equation |Q, |^. It is well known that for any integer value o = 1, 2, 3, . . ., the 
corresponding Lame potential (^) has a bound bands followed by a continuum band ||2|, ^. All band 
edge energies and wave functions are analytically known. 

At this point it is worth recalling that supersymmetric quantum mechanics (SUSYQM) has proved 
useful in discovering many, new, analytically solvable potentials on both the full as well as the half 
line Q. It is then natural to enquire if one can also use similar techniques to discover new solvable 
periodic potentials. In this paper, we demonstrate that this is indeed possible. 



Our work is inspired by several recent papers fs], |^, 0, ^ which discuss various general aspects 
of SUSYQM for periodic potentials. In particular, Dunne and Feinberg |^] defined and developed 
the concept of "self-isospectral" periodic potentials in detail. A one dimensional potential V-{x) of 
period L is said to be self-isospectral if its supersymmetric partner potential V+{x) is just the original 
potential upto a discrete transformation - a translation by any constant amount, a reflection, or both. 
A common example is translation by half a period, in which case the condition for self-isospectrality 
is 

V+ix)=V.ix-L/2) . (3) 

It is easily checked that if the superpotential W satisfies 

W{x) = -W(x - L/2) , (4) 

then condition (|3|) immediately follows. In this sense, any self-isospectral potential is rather uninter- 
esting, since application of the SUSYQM formalism Q to it just yields a discrete transformation and 
basically nothing new. We have recently pointed out ^J that the Lame potentials given in eq. (||) are 
not self-isospectral for a > 2, and hence SUSYQM generates new exactly solvable periodic problems. 
This point is further developed in detail in this paper. 

We expand our discussion to the band edges and wave functions of a much richer class of periodic 
potentials given by 

V{x) = pm sn^(x) -|- qm — jZT ' P — ^(^ + 1) ) <? = ^(& + 1) , (5) 

dn [x) 

where, like sn(2;), the Jacobi elliptic functions cn(x) and dn(x) also have a modulus parameter m 
which, for notational convenience, is not explicitly displayed. The potentials of eq. (|5|) are called 
associated Lame potentials, since the corresponding Schrodinger equation is called the associated 
Lame equation Q . More precisely, we often refer to the associated Lame potential of eq. (g) as the 
(p, q) potential and note that (p, 0) potentials are just the ordinary Lame potentials. Although some 
results for (p, q) potentials are available in scattered form in the mathematical literature, many of 
our results are new. In particular, we obtain band edge energies and wave functions for the special 
case p = q = a{a + 1) for a = 1, 2, 3, . . .. We study many (p, q) potentials and check whether they 
are self-isospectral by constructing and examining the supersymmetric partner potentials. In most 
cases, V-{x) is not self-isospectral, and consequently V+{x) is a new, exactly or quasi-exactly solvable 
periodic potential. 



The associated Lame potentials given by eq. (^) can also be re-written in the alternative form 

V{x) = pm sn (x) + qm sn {x + K[m)) , (6) 

since Hll] 



sn(x + K) = cn(x)/dn(x) , cn(x + K) = — Vl — m sn(x)/dn(x) , dn(x + K) = ^J\ — m/dn{x) . 

It is clear from (g) that potentials {p,q) and {q,p) have the same energy spectra with wave functions 
shifted by K{m). Therefore, it is sufficient to restrict our attention to p'> q. 

Before actually solving the Schrodinger equation for the associated Lame potential (g) , let us make a 
few general comments. Throughout this paper, we have chosen units with fi = 1, and taken the particle 
mass in the Schrodinger equation to be 1/2. Note that in the limit when the elliptic modulus parameter 
m = 0, the potential vanishes and one has a rigid rotator problem of period 2if (0) = vr, whose energy 
eigenvalues are aX, E = 0,1,4,9,... with all the nonzero values being two-fold degenerate. On the 
other hand, the limit m ^ 1 is much trickier since K{m) tends to infinity and the periodic nature 
of the potential is obscured. The Schrodinger equation for finding the eigenstates for an arbitrary 
periodic potential is called Hill's equation in the mathematics literature [y]. A general property of 
Hill's equation is the oscillation theorem which states that for a potential with period L, the band edge 
wave functions arranged in order of increasing energy Eq < Ei < E2 < E^ < E4 < E^ < Eq < ... are 
of period L, 2L, 2L, L, L, 2L, 2L, .... The corresponding number of wave function nodes in the interval 
L are 0, 1, 1, 2, 2, 3, 3, ... and the energy band gaps are given by Ai = E2 — Ei, A2 = -E4 — E^, A3 = 
Eq — E^, .... We shall see that the expected m = limit and the oscillation theorem are very useful 
in identifying if all band edge eigenstates have been properly determined or if some have been missed. 

The plan of the paper is as follows. In Sec. 2, we briefly review the basic ideas of SUSYQM. A 
detailed discussion of Lame potentials and their supersymmetric partners is given in Sec. 3. Solutions 
of the Schrodinger equation for the associated Lame potentials are presented in Sec. 4. Many key 
new results are summarized in Table 3. It is shown that the locus of quasi exactly solvable problems 
|lO| , prfl in the {p, q) plane are parabolas about the line p = q. Our solutions are valid for any real 
choice of the parameters a,b (recall p = a{a + l),q = b{b + 1)). Integer and half-integer values of 
a, b , including the very interesting special case a = b =integer, are treated in detail in Sec. 5. In 
most cases, the application of SUSYQM gives new solvable periodic potentials, many of which are 
illustrated in the figures. Finally, Sec. 6 contains some concluding remarks. 



2. Supersymmetric Quantum Mechanics Formalism: 

The supersymmetric partner potentials V±{x) are defined in terms of the superpotential W{x) by 

V±{x) = W^{x) ± W'(x) . (7) 

The corresponding Hamiltonians H± can be factorized as 

H_ = A+A, H+ = AA^ , (8) 

where 

A = ^ + W{x),A-^ = -^ + W{x), (9) 

so that the spectra of H± are nonnegative. It is also clear that on the full line, both H± cannot have 
zero energy modes since both t/j^ given by 

#)(x)=exp(±y'V(y)dy) , (10) 

cannot be simultaneously normalized. 

On the other hand, when the superpotential W{x) is periodic {W{x + L) = W{x)) then the 
potentials V-{x) and V+{x) are isospectral - their spectra match completely, including the zero modes, 
and SUSY is unbroken provided 

W{y)dy = 0. (11) 



It is worth noting that in this case both i/jq belong to the Hilbert space. Thus in this case even though 
SUSY is unbroken, the Witten index is zero Ig. The condition ([ll| ) is trivially satisfied in case W{x) 
is an odd function of x and throughout this paper we shall only consider superpotentials W which 
are odd function of x. Further, using the known eigenfunctions t/jn (x) of V-{x) one can immediately 
write down the corresponding un-normalized eigenfunctions ipn (x) of V+(x). In particular, from eq. 
( [ToD it follows that the ground state of V-^{x) is given by 



4~'\^) = ^^, (12) 



while the excited states t/jn (x) are obtained from t/jn (x) by using the relation 

i^i^Hx) = [^ + Wix)] V^(-)(x) , (n > 1) . (13) 



Thus by starting from an exactly solvable periodic potential V-{x), one gets another isospectral 
periodic potential V+(a;). As emphasized previously, if V-{x) is not self-isospectral, then V-^-{x) is a 
new solvable periodic potential! 

3. Lame Potentials (p,0) and Their Supersymmetric Partners: 

The supersymmetric quantum mechanics formalism of the previous section will now be applied to 
the Lame potentials ma{a + 1) sn^(2;, m). Analytic solutions are known for integer values of a 0, and 
the supersymmetric partner potentials can be readily computed. We first discuss the results for small 
integer values of a, and then present some eigenstate results for arbitrary integer values of a. 

A. Lame potentials with a = 1,2,3 : 

a=l: The a = 1 Lame potential VL. = 2m-sn^(x) — TTi is known to be self-isospectral |5| since its 
SUSY partner satisfies Vj^{x) =V-{x — K{m)). Both V"+(x) and V-{x) have one energy band ranging 
from energy to energy 1 — m, with a continuum starting at energy 1 [0]. Note that at m = one has 
energy eigenvalues at 0, 1 as expected for a rigid rotator and as ttt, — > 1, one gets V^{x) ^ 1 — 2sech x, 
the band width 1 — m vanishes as expected, and one has an energy level at E = 0. 

a=2: For the a = 2 case, the Lame potential (Q) has 2 bound bands and a continuum band. The 
energies and wave functions of the five band edges are well known |2|, |3|. The lowest energy band 
ranges from 2 + 2m — 26 to 1 + m, the second energy band ranges from 1 + 4m to 4 + ?n, and the 



continuum starts at energy 2 + 2m + 26, where 6 = Vl — m, + m,"^. The wave functions of all the band 
edges are given in Table 1. Note that in the interval 2K{m) corresponding to the period of the Lame 
potential, the number of nodes increases with energy. In order to use the SUSYQM formalism, we 
must shift the Lame potential by a constant to ensure that the ground state i.e. (the lower edge of 
the lowest band) has energy E = Q. As a result, the potential 

V_ {x) = -2 -2m + 26 + &msn^{x) (14) 

has its ground state energy at zero with a corresponding un-normalized wave function [Q] 

ipQ {x) = 1 + m, + 6 — 3m,sn {x) . (15) 

The corresponding superpotential is 

d , ,(-\, ■, 6m,sn(x)cn(3;)dn(3;) , , 

VF = -— log^^xH \ -, , 16 



and hence the partner potential V+(x) for the potential V-{x) given in eq. (jTJ) is 

y (^) _ y (^) I 72mhn^{x)cn^x)dn\x) 

V+[X)- ^^-W+ [i + ^ + ^_3^gn2(^)]2 • 1^0 

Although the SUSYQM formalism guarantees that the potentials V± are isospectral, they are not 
self-isospectral, since they do not satisfy eq. @ [^. Therefore, V+{x) as given by eq. (^) is a new 
periodic potential which is strictly isospectral to the potential (|l^ ) and hence it also has 2 bound 
bands and a continuum band. In Fig. 1 we have plotted the potentials V±{x) corresponding to a = 2 
for three different values of the parameter m. The values are m = 0.5, 0.8, 0.998. The difference in 



shape between V-{x) and V+(x) is manifest from the figures, especially for large m. Using eqs. (12) 



and (13) and the known eigenstates of V-{x), we can immediately compute all the band-edge Bloch 
wave functions for V+(x). In Table 1 we have given the energy eigenvalues and wave functions for the 
isospectral partner potentials V±(x). At m = one has energy eigenvalues 0,1,4 as expected for a 
rigid rotator. As rn, ^ 1, one gets V-{x) — > 4 — 6sech x, the band widths vanish as expected, and one 
has two energy levels at E = 0, 3, with a continuum above E = 4. 

a=3: For the a = 3 Lame potential, the ground state wave function is 

■ipQ (x) = dn(x)[2m + 5i + 1 — Bmsn (x)] , 

the corresponding superpotential is [^ 

?TT,sn(x)cn(x) [2?tt, + (5i + 11 — 15msn^(x)l , , 

vV = — - — - — - — -- (18) 

dn(x) [2m + (5i + 1 — 5msn^(x)] 

and the partner potentials V±{x) are Q 



V-{x) = -2 - 5m + 26i + 12m sn2(x) , 6i = Vl-m + Am^ , (19) 

and 

2m^sn^(x)cn^(x) [2m + 5i + 11 — 15m sn^(x)]^ 
14(x) = -l/_(x)+ ^-2^-^ [2m + 5, + l-5msn^{xW ' ^ ' 

Clearly, the potential V-{x) is not self-isospectral. In fact, y_(x) and V+{x) are distinctly different 

periodic potentials which have the same seven band edges corresponding to three bound bands and a 

continuum band ||2[. In Fig. 2 we have plotted the potentials V±(x) corresponding to a = 3 for several 

different values of the parameter m. The values of m are 0.5, 0.8, 0.998. It is clear from the figure that 

the potentials V+(x) and V^{x) have different shapes and are far from being self-isospectral. Using 

eqs. (|l2[) and (O) and the known eigenstates of V-{x), we can immediately compute all the 7 band 



edges corresponding to the known 3 bound bands and a continuum band |§, |3|. For example, the 
ground state -00 is given by 

1 1 

^° ^' 4"^(x) dn(x)[l + 2m + (5i-5msn2(x)] ■ ^ ^ 

The wave functions for the remaining six states are similarly written down by using eq. (^). These 
are shown in Table 2. The band edge energies for the a = 3 Lame potential (12,0) as a function of the 
elliptic modulus parameter m are plotted in Fig. 3. Note that at m = one has energy eigenvalues at 
0, 1, 4, 9 as expected for a rigid rotator and as m — > 1, one gets V- (x) ^ 9 — 12sech x, the band widths 
vanish as expected, and one has three energy levels at -B = 0, 5, 8 with a continuum above E = 9. 

B. Results for general integer values of a: 

The extension to higher values of a is straightforward. It is possible to make several general 
comments about the form of the band edge wave functions for the partner potentials V^{x). This is 
most conveniently done by separately discussing the cases of even and odd values of a. 

a= even integer: For a even, say a = 2N, it is known Q that there are A^ + 1 solutions of the 
form i<7v(sn2x), and N solutions each of the three forms 

snx cnx FAr_i(sn x) , snx dnx F]\[-i{sn x) , cnx dnx F/v-i(sn x) . 

Here Fr denotes a polynomial of degree r in its argument. The ground state i1^q{x) (which is the 
lower edge of the lowest band) is of the form F/v(sn2x). It is easily checked using eq. (|l3|) that the 
corresponding partner potential V+(x) has A^ solutions each of the four forms 

dnx Gj\f{sn'^x) snx GnIsu^x) cnx GN{sn^x) snx cnx dnxGAr_i(sn2x) 
-00" (x) ' ipQ{x) ' il^^ix) ' ^pQix) 

while the ground state is given by V'o'(^) = l/V'o'(^)- 

a= odd integer: For a odd, say a = 2N + 1, it is known [g| that the Lame potentials have A^ + 1 
solutions each of the three forms 

snx F/v(sn^x) , cnx F/v(sn^x) , dnx FAr(sn^x) 

and A^ solutions of the form 

snx cnx dnx F/v-i(sn x). 



The ground state ipQ (x) is of the form dnx Fj\f {sii^x). We can then easily deduce that the corresponding 
partner potentials Vj^{x) will have A^ + 1 solutions each of the two forms 

snx Gat+i (sn^x) cnx GAr_|_i(sn^x) 
i/joix) ' i^oix) 

and A'^ solutions each of the two forms 

dnx GAr+i(sn^x) snx cnx dnx G7v(sn^x) 



IpQ (x) ^0 (x) 

while as usual, the ground state is given by V'o~(x) = 1/iPq{x). 

In summary, for integral a, Lame potentials with a > 2 are not self isospectral. They have distinct 
supersymmetric partner potentials even though both potentials have the same (2a + 1) band edge 
eigenvalues. 

4. Associated Lame Potentials (p,q) and Their Supersymmetric Partners: 

In contrast to the Lame potentials discussed above, there seems to be no systematic treatment of 
associated Lame potentials in the literature. Therefore, we will first devote some time to discuss the 
properties of associated Lame potentials, show that they are quasi exactly solvable and then proceed 
to construct and study their isospectral supersymmetric partner potentials. 

A. Description of associated Lame potentials: 

As mentioned before, we will refer to the associated Lame potentials given by eq. (^) or the 
equivalent form eq. (g) as the {p, q) potential. The special cases p = 0, as well as q = 0, correspond 
to ordinary Lame potentials. 

In general, for any value of p and q, the associated Lame potentials have a period 2K{m) since 

sn(x + 2K) = — sn(x) , cn(x + 2K) = — cn(x) , dn(x + 2K) = dn(x) . 

However, for the special case p = q, eq. (|6|) shows that the period is K{m). From a physical viewpoint, 
if one thinks of a Lame potential {p, 0) as due to a one-dimensional regular array of atoms with spacing 
2K{m), and "strength" p, then the associated Lame potential {p,q) results from two alternating types 
of atoms spaced by K{m) with "strengths" p and q respectively. If the two types of atoms are identical 
[which makes p = q], one expects a potential of period K{m). 

Extrema (defined for this discussion as either local or global maxima and minima) of associated 
Lame potentials are easily found by setting dV{x)/dx = 0. This gives 

sn(x) cn(x) \p dn^(x) — g(l — m)] = . 



Extrema occur when (i) sn(x) = 0, that is x = 0, it 2K{m), it iK^m), . . . ; (ii) cn(j;) = 0, that is 
X = lb K{m), lb 3K{m), . . . ; (iii) dn (x) = (1 — m)q/p . At the points specified by (i) and (ii), one 
always has extrema and V{x) has values pm and qm. In addition, since dn (x) has a minimum value 
(1 — m)^ and a maximum value unity Q, condition (iii) yields additional extrema provided 

(1 — m) < (1 — m)q/p < 1 . 

For given fixed values of q and m, this condition has a solution provided p lies in the critical range 

g(l — m) < p < q/il — m) . 

Alternatively, for given fixed values of p and q with p > q, condition (iii) has a solution provided m is 
greater than the critical value 1 — q/p. 

The associated Lame potentials for q = 2,m = 0.5 and several values of p are plotted in Fig. 4(a). 
In the critical range of p values 1 < p < 4 , one expects additional extrema, and these are clearly 
seen in Fig. 4(a). In general the period is 2K(0.5) = 3.708 , but for p = q = 2, the period K(0.5) is 
evident. Note that as p increases, any given extremum changes character. For example, at x = 0, as 
p increases, one goes from a maximum to a local minimum to an absolute minimum. In Fig. 4(b) we 
have plotted associated Lame potentials for p = 4, g = 2 and several values of m. As expected from 
the above discussion, one always sees extrema at the points specified by conditions (i) and (ii), and 
additional extrema coming from condition (iii) are evident for m > 1/2. 

B. Solutions of the associated Lame equation - parabolas of solvability: 

The associated Lame equation is just the Schrodinger equation for the potential in eq. (|5|). 

d'^tl' r 9/ N cn^(x) _, , , , 

- Vt + b^ sn^ X + qm ^ ^) - g^ = . 22 

ax^ dn (x) 

On substituting 

^(x) = [dn(x)]-'' y{x) , (23) 

it is easily shown that y(x), satisfies the Hermite elliptic equation Q 

y"{x) + 2bm^^^^^^^^y'{x) + [X - {a + 1 - b){a + b)msn\x)]y{x) = , (24) 

dn(x) 

where 

p = a{a + l) , q = b{b + l) , E = X + mb"^ . (25) 

10 



On further substituting 

sn(x)=sint, y{x) = z{t) , (26) 

one obtains Ince's equation 

(1 -msin2t)/(t) + (26-l)msintcost z'{t) + [A- {a + 1 -h){a + h)m sir? t\z{t) =0, (27) 

which is a well known QES equation |^]. In particular, on substituting 

°° u"'R 
cos t = u , z(t) = w{u)w(u) = y^ -^ , (28) 

n=0 

it is easily shown that Rn satisfies a three-term recursion relation. In particular ifa + 6+1 = n ( 
n = l,2,3,...) then one obtains n QES solutions. Actually n QES solutions are also obtained in case 
b — a = —n{n = 1,2,3, ...) but since q is unchanged under 6 —> —6 — 1, no really new solutions are 
obtained in this case. The QES solutions for n = 1, 2, 3, 4, 5 are given in Table 3. In particular, for any 
given choice of p = a{a + 1), Table 3 lists the eigenstates of the associated Lame equation for various 
values of q. 

For q = a{a— 1), there is just one eigenstate with energy ma?' and wave function ^ = dn"'{x). Since 
the wave function has period 2K{m) and is nodeless, this is clearly the ground state wave function of 
the {a{a + l),a(a — 1)) potential for any real choice of the parameter a. The equations p = a{a + 1) 
and q = a{a — 1) are the parametric forms of the equation of the parabola {p — q)"^ = 2(p + g), which 
is plotted in Fig. 5 and labeled PI. For any point on the parabola, one knows the ground state wave 
function and energy Eq = ma?. The parabola PI includes the points (2,0) and (6,2). 

For g = (a — l)(a — 2), we see from Table 3 that two eigenstates with energies 1 + m{a — 1)^ and 
1 + 7Tia^ are known. Since they have period 4K(m,) and just one node in the interval L = 2K{m,), they 
must correspond to the first and second band edge energies Ei and E2 of the {a{a + 1), (a — l)(a — 2)) 
potential. Eliminating a from the equations p = a{a + 1) and q = (a — 1) (a — 2) gives the "parabola of 
solvability" {p — q)"^ = 8{p + q) — 12, which is plotted in Fig. 5 and labeled P2. This parabola includes 
the points (2,0) and (6,0) which correspond to Lame potentials. Similarly, the parabolas of solvability 
Pn (n = 0, 1, 2, ...) corresponding to q = {a — n + l)(a — n) in Table 3 are plotted, n eigenstates are 
known for any point on the parabola of solvability Pn. 

C. Supersymmetric partner potentials: 

It is easily checked from Table 3 that the solution corresponding to q = a{a — 1) as well as one 
of the q = (a — 2) (a — 3) solutions are nodeless and correspond to the ground state. Hence, for these 

11 



cases, one can obtain the superpotential and hence the partner potential V^ and enquire if V- is 
self-isospectral. For example, consider the case of p = a{a + l),q = a{a — 1) in which case W is given 
by 

^^_V;^^^^sn(x)cn(x)^ 
t/^oix) dn(x) 

so that the corresponding partner potentials are 

V- = [a — l)am — j—T + '>tT'(^{o- + Ijsn [x) — ma , 
dn (x) 



V"+ = a{a + l)m ^ ^ ^ + raia — l)asn (x) — ma . (30) 



cn^(x) 
dn (x) 

It is easily seen that these partner potentials satisfy eq. (y), are consequently self-isospectral and 

SUSY gives nothing new in this case. It is amusing to note that the superpotential W obtained here 

was in fact discussed in ref. [^ (see their eq. (32)). 

Let us now consider the SUSY partner potential computed from the ground state for the p = 

a{a + 1), g = (a — 2) (a — 3) case. It is given by (see Table 3) 

V'o(x) = [m(a-l) -l-5i+m(2a-l)sn2(x)](dn(x))''-2 ^ (g^) 

where 5i = yjl — tti + m^^{a — 1)-^. The corresponding superpotential W turns out to be 

m{a — 2)sn(x)cn(x) 2r?T,(2a — l)sn(x)cn(x)dn(x) 

~ dn(x) [m(l -a)-l-5i+ m{2a - l)sn2(x)] ' ^ ' 



Hence the corresponding partner potentials are 



cv?(x) 



V-{x) = ma{a + l)sn^(x) + m{a - 3)(a - 2) — ^^ - 2 - m{a^ - 2a + 2) + 26i , (33) 

dn (x) 

V+{x) = -V-{x) + 2W^{x) . (34) 

It is easily checked that these potentials are not self-isospectral since they do not satisfy the condition 
(g|). Thus one has discovered a whole class of new elliptic periodic potentials V+(x) as given by eq. 



(34) for which three states are analytically known no matter what a is. In particular, the energy 
eigenfunctions for V+ of these three states are easily obtained by using the corresponding energy 
eigenstates of V- as given in Table 3 and using eqs. (|l^ ) and (p^). 

5. Associated Lame Potentials with Special Values of p and q: 

We shall now discuss associated Lame potentials {a{a + 1), b(b + 1)), where a and b are either both 
positive integers or half- integers. In most cases, we show that although several band edge energies are 
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exactly known from Table 3, one usually does not know all the band edge energies, that is one has a 
quasi exactly solvable problem. However, in the special case of p = q{a = b = integer), we show that 
all the band edge eigenstates can be obtained and one has an exactly solvable periodic problem. 

A. a,b = integer, a 7^ b: 

First, let us note that the Lame potentials {a{a + 1),0) are in this category when a = integer 
and 6 = 0. For example, when a = 3, one has the (12,0) potential. We see from Fig. 5 that two 
parabolas of solvability pass through the point (12,0). From Table 3 it follows that 3 band edges of 
period 2K{m) are obtained from q = (a — 2) (a — 3) and 4 band edges of period 4:K{m) are obtained 
from q = (a — 3) (a — 4). Altogether, arranging in order of increasing nodes, one has 7 band edges 
with periods 2K, AK, 4K, 2K, 2K, AK, AK with 0, 1, 1, 2, 2, 3, 3 nodes respectively. There are no missing 
states, and as discussed in Sec. 3A, this gives three bound bands and a continuum band. 

As a second example with g 7^ 0, consider the (6,2) associated Lame potential, that \s p = Q,q = 2. 
In this case, taking a = 2, one can get five band edges from Table 3 - one solution of period 2K 
is obtained from q = a{a — 1), while the remaining four solutions of period AK are obtained from 
q = (a — 3) (a — 4). The eigenvalues and eigenfunctions are given in Table 4 along with the number 
of nodes in one period 2K. It is clear that there are two solutions of period 2K with 2 nodes in the 
interval 2K which have to be present but have not been obtained. This is also clear from the m = 
limit, since the energies from Table 5 are 0,1,1,9,9 and the states at 4,4 are missing. Thus, this is a 
QES problem. Fig. 6 illustrates the (6,2) associated Lame potential and its supersymmetric partner 
for three choices of m. The self-isospectral nature of the (6,2) potential is evident from Fig. 6 - it also 



follows from eqs. (30) with a = 2. The band edge energies for the (6,2) associated Lame potential as 
a function of the elliptic modulus parameter m is shown in Fig. 7. The two unobtained band edges of 
period 2K will have energies E = A at m = and E = 3 at m = 1. 

Let us now discuss the general associated Lame potential {a{a + l),b{b + 1)). Without any loss 
of generality let us assume that a > b. Using Table 3, we obtain (a — b) states of period 2K{AK) for 
q = [a — (a — 6)] [a — (a — 6 — 1)] for (a — b) odd (even), and {a + b + 1) states of period 4i^(2i^)for 
q = [a — {a + b + \)][a — {a + b)] for (a — b) odd (even). It can be established that some states are 
missing by looking at the node structure as well as the m = limit. Hence we again have a QES 
problem. 
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B. a = b = integer: 

Let us now discuss the special case oi p = q = a{a + 1), a = 1,2, ... . In this case the associated 
Lame potential (^) has period K, rather than 2K. It then follows from the oscillation theorem that 
with increasing energy, the band edges must have periods K, 2K, 2K, K,K, ... and in the m = limit 
the eigenvalues must go to i? = 0, 4, 16, 36, ... with all nonzero eigenvalues being doubly degenerate. It 
is easy to check from Table 5 that one case for which we already have exact results is when p = q = 2. 
In particular, consider the special case a = 1, for which V^{x) of eq. ( |34 ) takes the form 



V-{x) = 2msn^{x) + 2m ™„^^'^ - 2 - m + 2y/l - m . (35) 

dn (x) 

Using Table 5, we can calculate three energy eigenvalues and eigenfunctions of V- taking a = 1 in 
q = (a — 2) (a — 3). These are given in Table 5. Whereas the ground state is of period K, the next 
two states in Table 5 indeed have period 2K. Using a = 1 in eqs. (|3l|) to (|j), we find that the 
corresponding SUSY partner potential is 

8-v/l — mTTi^sn^ ( j;)cn^ (x) 
~\d^?{x)~+VT^mf 



V+{x) = 2-m-2Vl-m \^ ^ ^ — j^J==^ . (36) 



Are the potentials V±{x) self-isospectral? Using the relations 



sn(x + K{m)/2) = (1 + \/l - m)^^'^ 
cn(x + K{m)/2) = (1 + ^/l-mf^il - rr,\^l'^ 



m 



' \Jl — m sn(x) + cn(x)dn(x) 
dn (x) + \/l — m 

(1 + Vl ~ r?T,)^"cn(3;) — sn(x)dn(x) 



dn(x + K{m)/2) = (1 - m) 



1/4 



dn (x) + \/l — m 
1 + \/l ~ m)dn(a;) — msn(x)cn(x) 



dn (x) + ^/l — m 
a little algebra reveals that indeed V± are self-isospectral and satisfy eq. (^). 



(37) 

(38) 
(39) 



Are the higher members of the p = q family (i.e. p = q = 6,12, 20, ...) also self-isospectral? If our 
experience with the Lame case is any guide then we would doubt it. Indeed, we will now show that 
the (6,6) associated Lame potential is not self-isospectral. We get five band edges analytically from 
Table 3. In particular, take a = 2 and consider the case of q = (a — 4) (a — 5), for which we know 
two eigenstates as given in Table 3. In fact, in this case three more eigenstates can be analytically 
obtained but the corresponding eigenvalues and eigenfunctions have not been given in Table 3 since 
the energy eigenvalues are solutions of a cubic equation whose exact solution for arbitrary a can not 
be written in a compact form. However, for a = 2, we are able to solve the cubic equation and obtain 
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the three eigenvalues in a closed simple form. In particular consider an ansatz of the form 

y = A + Bsn^x + Dsn^x . (40) 

On substituting this ansatz in eq. (^) it is easy to show that the energy eigenvalue A(= E — m{a — 4:)'^) 
must obey the cubic equation 

A^ + [28m - 20 - 12am] A^ + [64 - 304m + 160ma + 32m^ (a -2) (a -3)] A -64m(2a- 3) (2- 2m + ma) =0. 

(41) 
The solution of this equation is in general quite lengthy but in the special case of a = 2 this cubic 
equation is easily solved yielding three eigenvalues in a compact form. On combining them with the 
two levels given in Table 3, we obtain the eigenvalues and eigenfunctions of all the five band edges 
for the case p = q = 6. These are given in Table 6. We have also verified that these five eigenstates 
in ascending order of energy indeed have periods K, 2K, 2K, K, K respectively and taht the energy 
eigenvalues have expected limits at ?ti = 0. In particular the associated Lame potential V-{x) is 



cn^(2;) 
dn^(x) 



V-{x) = 6msn2(x) + 6m ^ ^ [ - 8 - 2m + 2^8 , (42) 



whose ground state energy is zero while the corresponding eigenfunction iJ)q is 
1 - (4 - m - 5^)s-n^{x) + (4 - 2m - 5s)s-a^{x) 



IpQ (X) 



dn^(x) 

Hence the corresponding superpotential is 

, , — 2msn(x)cn(x) 2sn(2;)cn(x) 
W{x) = \ \ ^ ^ + ^ ^ ^ ^ 



58 = VlG - 16m + m2 . (43) 



dn(x) dn(x)V'o {x 

and the partner potential V+{x) which is isospectral to V-{x) is 



(4 - m - Js) - 2(4 - 2m - 5s)sn^ix) 



(44) 



V+{x) = -V-{x) + 2W'^{x) . (45) 



It is not difficult to see that the W as given by eq. ( [441) does not satisfy the self-isospectral condition 
@ and hence unlike the p = q = 2 case, the p = q = 6 potential is not self-isospectral. In Fig. 8, 
we have plotted the potentials V±{x) corresponding to p = q = 6 for several different values of the 
parameter m. The figures confirm that the potentials are far from being self-isospectral. Thus we have 



obtained a new exactly solvable periodic potential (45) which has two bound bands and a continuum 



band, with five band edges and the corresponding eigenfunctions being exactly known using Table 6 
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and eqs. (p^ and (|13|). In Fig. 9, we plot the band edge energies for the (6,6) potential as a function 
of the elliptic modulus parameter m. 

It is also clear from here that even the higher associated Lame potentials with p = q = 12,20, ... 
which have 7,9,... band edges are also exactly solvable in principle and none of them will be self- 
isospectral, so that in each case one obtains a new exactly solvable periodic potential. In particular, 
for p = q = n{n + 1) there will be {2n + 1) band edges in both V±{x) whose energy eigenvalues can be 
obtained from Table 3 when q has the form [n — 2n][n — (2n + 1)]. Out of the (2n + 1) band edges in 
V-{x), (n + 1) solutions (including the ground state) have the form dn"a:^ while n solutions have the 
form F„_i(sn^x)^^^^^. On the other hand, as far as the (2n + 1) solutions of the partner potential 
V+ are concerned, there are n states each of the two forms 

snxcnxG„(sn^x) G„+i(sn^x) 



while the ground state (i.e. the lower edge of the lowest band) is given by iPq{x) = 1/iPq{x). 

C. a,b = half-integer: 

Let us now specialize to the case when both a, b are half integral with a > b. As an illustration, 
let us first consider the case of a = 3/2, b = 1/2 so that p = 15/4, q = 3/4. In this case, the oscillation 
theorem requires band edges with periods 2i^, 4i^, 4i^, 2i^, 2i^, .... Using Table 3 and Fig. 5, we see 
one gets three eigenstates when q = (a — 2) (a — 3) with a = 3/2, all with period 2K. The ground 
state is at £^0 = ^ while there are two degenerate levels at E^ = E4 = 4 -\- ^ . To understand this 
degeneracy better, let us go along the parabola of solvability P2 given hy q = (a — 2) (a — 3). The 



band gap is given by A2 = | —2 + m + 2y^l — m + 'm?{a — 1)^1 and is plotted in Fig. 10. It vanishes at 
a = 3/2 (15/4, 3/4) potential, and has the correct values A2 = 2\/l — m + m? — 2 + m for a = 2 and 



A2 = 2\/l — m + ^m? — 2 + 171 for a = 3 which correspond to the (6,0) and (12,0) Lame potentials. 
The vanishing of A2 at a = 3/2 occurs because the eigenfunctions corresponding to E^ and E4 cross 
over as one goes along the parabola P2. 

These arguments are easily generalised in case p = (n + 1/2) (n + 3/2), q = {k + 1/2) (/c + 3/2) with 
n > k. The energy eigenvalues of (n — k) states can be obtained by using Table 3 in case q is of the 
form q = [n + 1/2 — (n — k)][n + 1/2 — (n — k — 1)] and the corresponding eigenstates have period 
2K{4K) depending on whether (n — k) is odd (even). On the other hand, the energy of (n + /c + 2) 
states is obtained when q is of the form q = [n + 1/2 — (n + k + 2)][n + 1/2 — (n + /c + 1)] and these 
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states have the same period 2K(4K) as the n — k states when n — k is odd (even). It turns out that 

the n — k solutions are in fact common in both and so we only obtain the energy of the n + k + 2 band 

edges and all of them have the same period 2K{4:K) depending on if n — A; is odd (even) so that it is 

only a QES problem and not an exactly solvable problem as one is unable to obtain a single eigenstate 

with period 4K{2K) in case n — A; is odd (even). 

We would like to point out some of the pecularities of the spectrum in these cases. For example, 

in case {p,q) = (35/4,3/4), (63/4, 3/4), (99/4,3/4)... then one finds that 3,4,5,... QES energy levels 

of period 4K,2K,4K, ... respectively are analytically known of which the one at the highest energy 

is doubly degenerate. As an illustration, in Table 7 we have given the 4 QES energy eigenstates all 

of period 2K for the (63/4, 3/4) potential. The interesting point about this case is that the partner 

potentials V±{x) are not self-isospectral and hence one has discovered a new QES potential where 4 

band edges of period 2K and the corresponding eigenfunctions are explicitly known. Of these, the one 

at ii^ = 16 + ^ is doubly degenerate, again due to crossover of energy levels. Using the ground state 

wave function, the superpotential is computed to be 

3msn(x)cn(x) 24msn(x)cn(x)dn(x) 

W = — (46) 

2 dn(x) [12msn2(2;) - 2 - 5m - V4 - 4m + 25m2] 

Using eqs. (^) and (|l3|) the eigenstates of the SUSY partner potential V+ are then determined. 

6. Comments and Conclusions: 

In this paper, we have discussed solutions of the type given in Table 3, which correspond to the 
parabolas of solvability shown in Fig. 5. Lame potentials {p, 0) with p = a(a + 1) and integer a, always 
have two parabolas of solvability passing through - one parabola gives all states of period 2K and the 
other gives all states of period AK. This provides a deeper understanding of why such Lame potentials 
are fully solvable [|l^]. Similarly, we have obtained eigenstates for a large class of associated Lame 
potentials {p,q). Further, using the formalism of supersymmetric quantum mechanics, we have been 
able to discover many new exactly solvable and quasi exactly solvable periodic potentials involving 
Jacobi elliptic functions. This is a very substantial improvement over the currently known small 
number of exactly solvable periodic problems. 
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Table Captions 

Table 1: The eigenvalues and eigenfunctions for the 5 band edges corresponding to the a = 2 
Lame potential V- which gives {p,q) = (6,0) and its SUSY partner V+. Here B = 1 + m + 6 and 



6 = Vl — m + rrfi. The potentials V± have period L = 2K{m) and their analytic forms are given by 
eqs. ( p^ and (|l7|) respectively. The periods of various eigenfunctions and the number of nodes in the 
interval L are tabulated. 

Table 2: The eigenvalues and eigenfunctions for the 7 band edges corresponding to the a = 3 Lame 
potential F__ which gives {p,q) = (12,0) and its SUSY partner Vj^. Here 5i = Vl — m + Am?; 82 = 
\/4 — m + m? ; (^3 = \/4 — 7m + Am? . The potentials V± have period L = 2K{m) and their analytic 
forms are given by eqs. (19) and ( [20| ) respectively. The periods of various eigenfunctions and the 



number of nodes in the interval L are tabulated. 

Table 3: Eigenvalues and eigenfunctions for various associated Lame potentials (p, q) with p = 
a{a + 1) and q = {a — n + l)(a — n) for n = 1,2, 3, .... The periods of various eigenfunctions and the 
number of nodes in the interval 2K{m) are tabulated. Here 5^ = yjl — m + m^{a — 1)'^ ; ^5 = 
\/4, — 7m + 2ma + m^{a — 2)^ ; 5^ = y/A — tti — 2m,a + m^^{a — 1)^ ; 5^ = ^^9 — 9m, + m!^{a — 2)'^ . 



Table 4: The five eigenvalues and eigenfunctions for the self-isospectral associated Lame potential 
corresponding to a = 2, 6 = 1 which gives (p, g) = (6,2). The potential is V-{x) = Qmsv?{x) + 
2m^^-4r\ — 4?Ti, and has period 2K{rn). The number of nodes in the interval 2K{m) is tabulated. 

Table 5: The three eigenvalues and eigenfunctions for the associated Lame potential corresponding 
to a = 6 = 1 which gives {p,q) = (2,2). The potential has period K{m,) and the number of nodes in 
the interval K{m) is tabulated. 

Table 6: The five eigenvalues and eigenfunctions for the associated Lame potential corresponding 
to a = b = 2 which gives {p,q) = (6,6). Here Ss = \/l6 — 16m + tti^ . The number of nodes in one 
period K{m) of the potential is tabulated. 

Table 7: Energy eigenvalues and eigenfunctions for the associated Lame potential corresponding 
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to a = 7/2,6 = 1/2 which gives {p,q) = (63/4,3/4). Here 5g = \/4 - 4m + 25m2 ; V-{x) = 

^msn^^x) + 4m j^^^ — 2 j^ + (^g. The last column gives the number of eigenfunction nodes in one 

period 2K{m) of the potential. 



Figure Captions 

Fig. 1: The (6,0) Lame potential V-{x) corresponding to a = 2 [thick line] as given by eq. (14) 
and its supersymmetric partner potential V+(x) [thin line] as given by eq. (|l7|) for three choices of m 
(a) 0.5 (b) 0.8 (c) 0. 



Fig. 2: The (12,0) Lame potential V-{x) corresponding to a = 3 [thick line] as given by eq. ([l9|) 
and its supersymmetric partner potential V+(x) [thin line] as given by eq. (pO) for three choices of m 
(a) 0.5 (b) 0.8 (c) 0.998. 

Fig. 3: Band edge energies for the (12,0) Lame potential corresponding to a = 3 as a function of 
the elliptic modulus parameter m. This figure is drawn using the eigenvalues given in Table 2. 

Fig. 4: (a) Plots of the (p, (?) associated Lame potentials for q = 2,m = 0.5 and several values of 
p. (b) Plots of the {p, q) associated Lame potentials for p = 4, g = 2 and several values of m. 

Fig. 5: Parabolas of solvability. This figure illustrates all associated Lame potentials (p, q) which 
are quasi solvable. Each parabola corresponds to a choice of q in Table 3. Parabola Pn is for 
q = (a — n + l)(a — n) for n = 1,2, 3, ..., and one knows n eigenstates for any point on it from Table 
3. 

Fig. 6: The (6,2) associated Lame potential V^{x) [thick line] and its supersymmetric partner 
potential V+(x) [thin line] for three choices of m (a) 0.5 (b) 0.9 (c) 0.998. 

Fig. 7: Band edge energies for the associated Lame potential (6,2) as a function of the elliptic 
modulus parameter m. This figure corresponds to Table 4. 
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Fig. 8: The (6,6) associated Lame potential V-{x) [thick line] as given by eq. ( |4^ ) and its 
supersymmetric partner potential V+(x) [thin line] as given by eq. ( |45D for three choices of m (a) 0.5 
(b) 0.9 (c) 0.998. 

Fig. 9: Band edge energies for the associated Lame potential (6,6) as a function of the elliptic 
modulus parameter m. This figure corresponds to Table 6. 

Fig. 10: Energy gap A2 = IS4— -E3I as one moves along the parabola of solvability P2 corresponding 
to q = {a — 2)(a — 3) and p = a{a + 1). 
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Table 1: Energy Eigenstates for V± Corresponding to a = 2 



E 



V.(-) 



[B - 3m sn2(x)]V'^+) 



Period Nodes 



m + 1 + 6 — Smsn^ (x) 
26 — 1 — m cn(a:)dn(x) 

26 — 1 + 2m sn(x)dn(x) 

26 + 2 — m, sn(a;)cn(a;) 



4(5 



m,+ 1 — 6 — 3msn^(x) 



1 



sn(a:)cn(a:)dn(a;) 



2K 





sn(x) [6m - {m + 1)B + msii^ (x) {2B - 3 - 3m)] 4K 1 

cn{x)[B + msn2(x)(3 - 2B)] 4K 1 

dn(x) [B + sn^ (x) (3m - 2B)] 2K 2 



2K 



Table 2: Energy Eigenstates for V± Corresponding to a = 3 



E 



i^(-) 



^-^M 



Period Nodes 







dn(x)[l + 2m, + 6i — 5r7T,sn^(x)] 
cn(x)[2 + m + 62 — 5msn^(x)] 



3 - 3m + 2(5i - 262 

3 + 2(5i - 2(53 sn(x)[2 + 2m + (53 - 5msn2(x)] 



2 - m + 2(5i 

4(5i 

3 - 3m + 2(5i + 2(52 



sn(x)cn(x)dn(x) 

dn(x)[l + 2m — 61 — 5msn^(x)] 

cn(x)[2 + m — 62 — 5m,sn^(x)] 



3 + 2(5i + 2(53 sn(x)[2 + 2m - (53 - 5msn2(x)] 



1 

10m(l — 771 + (52 — (5i)sn(x)cn^(x)dn (x) 

^ ' cn(a;)dn(a;) 

10m(l + (53 — 5i)cn(x)sn^(x)dn (x) 

-(i - 2msn {X)) s„(^)dn(a;) 

dn^(x)[l + 2m + (5i + (m - 2 - 2(5i)sn2(x)] 

sn(x)cn(x)dn (x) 

10r7T,(l — 771 — (52 — 5i)sn(x)cn^(x)dn (x) 

^ '^ cn(a;)dn(a;) 

10m(l — (53 — (5i)cn(x)sn^(x)dn^(x) 

-U-2msn (^XJJ sn(^)dn(a:) 



2K 
AK 

AK 



AK 



2K 


2 


2K 


2 


AK 


3 
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Table 3: Some Eigenstates for Various Associated Lame Potentials 



E 



dn "(x)^ 



Period Nodes 



a{a — 1) ma 

(a-l)(a-2) l + m(a-l)2 

(a-l)(a-2) l + ma^ 

(a-2)(a-3) 2 + m(a2 - 2a + 2) ± 2^4 

(a-2)(a-3) A + m{a-lf 

(a-3)(a-4) 5 + m(a2 - 4a + 5) ± 2(55 

(a-3)(a-4) 5 + ^(a^ _ 2a + 2) ± 2(^6 

(a-4)(a-5) l{) + m{a^ - Aa + b) ± 25j 



1 

cn(a::) 
dn(a::) 
sn(a;) 
dn(a;) 
[m(2a— l)sn^(a;) — 1+m— ma±54] 

sn(a::)cn(a::) 

di?(x) 

cn(a;)[m(2a— l)sn^(a;)— 2+2m— mai^s] 

di?(^0 
sn(a::)[m(2a— l)sn^(a;)— 2+m— ma±<56] 

sn(3;)cn(x)[m(2a— l)sn^(a;)— 3+2m— ma±<57] 
dn-*(x) 



2i^ 





4/i: 


1 


4/i: 


1 


2K 


2,0 


2K 


2 


AK 


3,1 


AK 


3,1 


2K 


4,2 



Table 4: Energy Eigenstates for the (6, 2) Potential 



E 



ilj^- 



Period Nodes 







5 - 3m - 2V4 - 3m 
5 — 2m — 2\/A — 5m + tti^ 



5 — 2m + 2\/4 — 5m + tti^ 
5 - 3m + 2V4 - 3m 



dn (x) 

!^I^l3msn^{x) - 2 - V4 - 3m] 



dn(a;) 

dn(x) 
sn(x) 



||[3msn^(x) -2 



TTT. — \/4 — 5m + 777.^1 



dnfir) [3m.Sn^(x) — 2 — TTT. + \/4 — 5m + tti^] 
gg[3msn^(x)- 2 + ^/4^3^ 



jm 



2i(: 





4i(: 


1 


AK 


1 


AK 


3 


AK 


3 



Table 5: Energy Eigenstates for the (2, 2) Potential 





E 




dn(a;)-(/;*-^'' 


Period 


Nodes 











K 
2K 
2K 






dn^(x) + \/l -"I 







4Vl-m 


dn (x) — \Jl — 771 
sn(x)cn(x) 


1 


2- 


-m + 2Vl - 


- 771 


1 



23 



Table 6: Energy Eigenstates for the (6, 6) Potential 



E dn {x)ip^ ' Period Nodes 

1 - (4 - m - (58)sn^(x) + (4 - 2m - (58)sn''(a;) 
-4 + 2m + 2(58 l-2sn'^{x) +msn^{x) 

2 — m — 6-v/l — m + 25s sn(x)cn(a;)[l — (1 — \/l — r7T,)sn^(a;)] 

2 — 771 + 6-v/l — m + 2(58 sn(a:)cn(x)[l — (1 + \/l — r7T,)sn^(a;)] 

4(58 1 - (4 - m + (58)sn2(a;) + (4 - 2m + 58)sn''(a;) 
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Table 7: Energy Eigenstates for the (63/4, 3/4) Potential 
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2 - m + (59 


sn(x)cn(a;)dn^(x) 
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[12msn2(x) - 2 - 5m + (59]dn2(x) 


2i^ 


2 


14 - 7m + (59 


sn(a;)cn(x)[l — 2st?{x)\ 


2K 


4 


14 - 7m + (59 


[1 - 8sn2(x)cn2(x)] 


2K 


4 
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